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ABSTRACT 

PG 1553+113 is a known BL Lac object, newly detected in the GeV-TeV energy range by H.E.S.S 
and MAGIC. The redshift of this source is unknown and a lower limit oi z > 0.09 was recently 
estimated. The very high energy (VHE) spectrum of PG 1553+113 is attenuated due to the absorption 
by the low energy photon field of the extragalactic background light (EBL). Here we correct the 
combined H.E.S.S and MAGIC spectrum of PG 1553+113 for this absorption assuming a minimum 
density of the evolving EBL. We use an argument that the intrinsic photon index cannot be harder 
than r = 1.5 and derive an upper limit on the redshift of z < 0.69. Moreover, we find that a redshift 
above z = 0.42 implies a possible break of the intrinsic spectrum at about 200 GeV. Assuming that 
such a break is absent, we derive a much stronger upper limit of z < 0.42. Alternatively, this break 
might be attributed to an additional emission component in the jet of PG 1553+113. This would be 
the first evidence for a second component is detected in the VHE spectrum of a blazar. 
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1. INTRODUCTION 

PG 1553+113 was discovered in the Palomar-Gree n 
survey of UV-excess stellar objects (|Green et al.lll986l ). 
It is classified as a BL La c object based on its featureless 
spectrum (Mil ler fc Gre en 1983; Falomo & Treves 1990) 
and its significant optical variability (Mil ler et al.lll988f ). 
PG 1553+113 is well studied from the radio to the X-ray 
regime. Based on its broad-band spectral energy distri- 
bution (SED), PG 1553+11 3 is now classified as a high- 
frequency peaked BL Lac (jGiommi" "ei~an[T995h . simi- 
lar to most of the AGNs detected at very high energies 
(VHE). 

The redshift of PG 1553+113 is essentially un- 
known. It w as initially determined to be z=0.36 
(iMiller fc Green 1983) b ut la ter this claim was with- 
drawn (|Falomo fc Treveg|1990[ ). To date no emission or 
absorption lines have been measured despite several ob- 
servation campaigns with optical instruments. No host 
galaxy was resolved in Hubble Space Telescope (HST) im- 
ages of PG 1553+ 113 taken during th e HST survey of 110 
BL Lac objects (|Scarpa et al.ll2000l ). The HST results 
were used to set a lower limit of z> 0.78 (jSbarufatti et al.l 
I2005D for PG 1553+113. A more recent publication 
claims a lower limit on the redshift of z>0.09 using the 
ESO-VLT optical spectroscopy survey of 42 BL Lacer - 
tae objects of unknown redshift (jSbarufatti et al.ll2006D . 
The possibility of a large redshift is of critical impor- 
tance to VHE observations due to the absorption of VHE 
photons (IStecker et all Il992t iHauser fc Dwekl [200 L ) by 
pair-production on the extragalactic background light 
(EBL). This absorption, which is energy dependent and 
increases strongly with redshift, distorts the VHE en- 
ergy spectra observed from distant objects. It is be- 
lieved that VHE photons above E > 100 GeV from ob- 
jects at distances z > 1 are completely absorbed. Note 
that the most distant objects detected so far at VHE 
are lES 1218+304 (z = 0.182, Albert et al. (2006a)) and 
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lES 1101-232 (z = 0.186. lAharonian et all (|2006af )). 

Recently, VHE 7-ray emission from PG 1553+ 113 was 
measured by H.E.S.S (lAh aroman et al.l l2006b) and by 
MAGIC (| Albert et a"ni2006H) . Both collaborations re- 
ported a soft energy spectrum with a differential pho- 
ton index of F = 4.0 ± 0.6 and F = 4.2 ± 0.3 respec- 
tively. The VHE data from the two measurements were 
used independently to derive an upper limit on the red- 
shift of the source of z < 0.74. This limit is based 
on the assumption of a minimum possible level of the 
EBL and the maximum hardness of the intrinsic VHE 
spectrum. Although the intrinsic VHE 7-ray spectra 
of the AGNs are not well known, it can be assumed 
that the intrinsic pho ton index of the source s is not 
harder than Pint = 1-5 (jAharonian et al.ll2006ah . How- 
ever, there are emission scenarios where the maximum 
possible VHE pho t on in dex is even harder (Pint =2/3, 
ikatarzvhski et al.l (|2006D V 

Here we use the combined H.E.S.S. and MAGIC spec- 
trum of PG 1553+113 to derive upper limits on its red- 
shift, assuming the two limits on the hardness of the 
intrinsic photon index (Fint = l-5 and Pint = 2/3). 

In addition, we present an alternative method to es- 
timate an upper limit on the redshift of PG 1553+113 
assuming that there is no break in the intrinsic VHE 
spectrum of the source. 



2. MINIMUM EBL 

The most likely reaction channel in the interaction of 
VHE 7-rays with the low energy photons of the EBL is 
pair production 7vhe + 7ebl e+e~. This reaction 
has its largest cross section when the center of mass en- 
ergy is roughly 2 times the threshold energy of 2meC^. 
The relevant EBL wavelength range for the absorption 
of VHE 7-rays spans from UV light (0.1 /im) to far in- 
frared (few 100 /im). This light predominantly consists 
of redshifted star light of all epochs and reemission of a 
part of this light by dust in galaxies. The intrinsic (de- 
absorbed) VHE photon spectrum, dN/dEi, of a blazar 



2 



E 
5 




I 

T 
T 

+ 
A 



T 



T 

A 
O 
O 
X 

▲ 

□ 
□ 



X (\i.m) 



Bernstein et al. 2002, 2005 
Brown et al. 2000 (HST/STIS) 
Cambresy et al. 2001 (piRBE/2MASS) 
Dole et al. 2006 (SPITZER) 
Dwek & Arendt 1998 (DIRBE) 
Dwek & Arendt UL 1998 (DIRBE) 
Edelstein et al. 2000 (Shuttie/UVX) 
Elbaz et al. 2002 (ISO) 
Fazio et al. 2004 (SPITZER) 
Finkbeiner et al. 2000 (DIRBE) 
Frayer et al. 2006 (SPITZER) 
Gorjian et al. 2000 (DIRBE) 
Hauser et al. 1 998 (DIRBE/FIRAS) 
Hauser et al. UL 1998 (DIRBE/FIRAS) 
Kashlinsky et al. 1996 
Kashlinksy & Odenwald 2000 
Lagache et al. 2000 (DIRBE) 
Lagache et al. UL 2000 (DIRBE) 
Madau & Pozzetti 2000 
Martin et al. 1991 (Shuttie/UVX) 
Matsumoto et al. 2005 (IRTS) 
Metcalfe et al. 2003 (ISO) 
Papovlch et al. 2004 (SPITZER) 
Toller 1983/Lelnert 1998 
Wright & Reese 2000 (DIRBE) 
Aharonlan et al. UL 2006 (HESS) 
Knelske et al. (low limit) 2004 



Fig. 1. — Energy density of the cxtragalactic background light (EBL). Direct measurements, galaxy counts, low and upper limits are 
shown by different symbols as described in the legend. The dashed black line represents the upper limit set by H.E.S.S. (lAharonian et all 
I200 6al). The black solid curve is the minimum EBL spectrum taken from .Kneiskc et al.. (.2004 ) for z=0. The model takes evolution of the 
EBL into account. 



located at redshift z can be determined using: 

dN/dEi = dN/dEobs x exp[r-y^(i?, z)], 

where dN/dEobs is the observed spectrum and Try^{E, z) 
is the optical depth. While the cross-section of the pair 
production is well known the spectral energy density 
of the EBL is not, which makes the calculation of the 
optical depth uncertain. Direct measurements of the 
EBL are difficult because of the strong foreground emis- 
sion, mainly consisting of reflected sunlight and ther- 
mal emission from zodiacal dust particles. Indirect up- 
per limits on the EBL density were obtained using fluc- 
tuation analyses of the detected back ground radiation 
(jKashhnskv et all 119961 : iKashlinskv fc Odenwaldii2000l ). 
Recently, an upper limit on the EBL density in the range 
from 0.8 to 3 /im was det ermined using measured VHE 
spectra of distant AGNs (|Aharonian "etaI1l2006aD . On 
the other hand, a robust lower limit o n the EBL is set by 
galaxy counts in different wavebands (iMadau k, Pozzetti 
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galaxy counts m d iiterent waveoand s (liviaaau &l rozzett 
20001 iFazio et al.ll2004l: lElbaz et al.ll2002t iMetcalfe et al 



2003tlPapovich et al.ll2004D . A summary of current direct 
and indirect measurements of EBL is given in Fig. [T] 

In the present study, we want to use the lowest possi- 
ble realistic level of the EBL in order to derive a mini- 
mum redshift, above which a break in the VHE spectrum 
of PG 1553-1-113 becomes evident. Several EBL models 
have been developed, taking into account the evolution 
of the EB L, which is rel a ted to star and galaxy evol u- 
tion (e.g. [Kneiske et a j] (|2004); Pri mack et all pOOS); 
IStecker et al.l (|2006[ )). Though different in approach, 
the models agree in their predictions within a factor of 
2. For this s t udy, w e use the lower limit model from 
iKneiske et alj ()2004[ ). which is just at the level of the 
direct lower limits set by the galaxy counts (see Fig.[T]). 

3. COMBINING THE SPECTRA OF PG 1553+113 
FROM H.E.S.S. AND MAGIC 

The differential VHE 7-ray spectra of PG 1553+113 
pubhshed by H.E.S.S. and MAGIC are shown in Fig. H 
The photon fluxes are multiplied by E"^ , which is equiva- 
lent to a vF{v) representation. The H.E.S.S. data are 
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Fig. 2. — Differential measured energy spectrum of PG 1553+113 
multiplied by E^ to represent energy density. The spectrum mea- 
sured by H.E.S.S. is shown by grey points, whereas the MAGIC 
spectrum is shown by black boxes. The last point in the MAGIC 
spectrum is a 95% upper limit. 



shown as grey circles and the MAGIC data by black 
boxes. For H.E.S.S., the data were taken in spring- 
summer 2005. In case of MAGIC, it is a combined data 
set taken between spring 2005 and winter 2006. The 
last spectral point is a 95% upper limit and will not be 
used in the further analysis. AGNs are known to be vari- 
able sources in flux so that, in general, it is not trivial 
to combine non-simultaneous data sets. In the present 
case, however, the agreement between the spectra mea- 
sured by H.E.S.S. and MAGIC is very good. As a test, 
we fitted a power law 



dN/dE = NqE 



-r 



(1) 



to the observed spectra in the overlapping energy re- 
gion between 150 and 600 GeV by fixing the slope 
(observed photon index F) to the value measured by 
H.E.S.S.: r = 4. The resulting normalization factors 
are iVo = (2.01 ± 0.34) x lO^^^ phot/(TeV cm^ s)-i and 
No = (2.24 + 0.59) x 10"" phot/(TeVcm2 s)"^ for the 
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Fig. 3. — Constraint on the redshift of PG 1553+113. Hollow 
points: measured combined (normalized) differential energy spec- 
trum of PG 1553+113 using MAGIC (squares) and H.E.S.S. (cir- 
cles) data from 2005 and 2006. Filled points: intrinsic spectrum 
of PG 1553+113, using minimum possible density of the evolving 
EBL and the redshift of z = 0.42. Black dotted line: power law fit 
to the intrinsic spectrum; the fitted photon index is listed in the 
inlay. Black solid line: curved power law fit. The probability of 
the likelihood ratio test is 95.06%, meaning that the fit by a curved 
power law is to be preferred over the fit by a power law. 



H.E.S.S. and MAGIC data respectively. These numbers 
differ by 10%, which is small compared to the individ- 
ual statistical errors. In order to increase the statistical 
power of the tests described below, we use the combined 
spectrum of MAGIC and H.E.S.S. data on PG 1553+113. 
We normalize the spectra by multiplying the H.E.S.S. 
fluxes by a factor 1.1 to avoid a possible bias in our re- 
sults. In order to show the effect of the normalization 
we also perform the same study on a combined spectrum 
without this normalization. 

4. DETERMINATION OF UPPER LIMITS ON THE 
REDSHIFT OF PG 1553-^113 

Two different methods are used. The first one assumes 
that the spectral index of the intrinsic differential photon 
spectrum is not higher than 1.5 (or, alternatively, 2/3). 
The second one assumes that the intrinsic photon spec- 
trum has no break in the measured energy region from 
80 to 600 GeV. 

MAXIMUM INTRINSIC PHOTON INDEX Tint 

A fit of Eq. [2] is performed to the intrinsic spectrum 
of PG 1553-1-113. The fitted index Tint and its error ur 
are combined to Fmax, where Fmax — Tint + 2 x ap corre- 
sponding to a 95% confidence level. We note that despite 
the fact that for z > 0.42 (see below) a curved power law 
gives a significantly better fit to the intrinsic spectrum 
than a simple power law, the value of the latter fit is 
still acceptable. A redshift z is considered to be unre- 
alistic if Fmax < 1.5 or, in case of the extreme scenario, 

Fmax ^ 2/3. 

PRESENCE OF A BREAK IN THE INTRINSIC 
PHOTON SPECTRUM 

The second method is based on the indication that 
for larger redshifts (z > 0.3) the intrinsic spectrum of 
PG 1553-1-113 exhibits a break at about 200 GeV with 



the intrinsic energy spectrum rising after the break (see 
Fig. [3|). Such a break in the intrinsic spectrum of 
PG1553H-113 cannot be excluded a priori. It could at- 
tributed to a second population of 7-ray emitting par- 
ticles. However, in none of the measured VHE 7-ray 
spectra of extragalactic sources such a component has 
been found. Thus, if a spectrum shows a break, either it 
is the first time a second emitting component is found in 
a VHE 7-ray spectrum or a lower redshift value z has to 
be assumed. 

To test the presence of a second component (or 
the presence of a break) in the VHE spectrum of 
PG 1553-1-113 we performed a likelihood ratio test 
(jEadie et al.lll998l ) on the intrinsic spectrum. Two hy- 
potheses are tested. Hypothesis A is a simple power law 
(2 free parameters): 



dN/dE = NqE- 



(2) 



Hypothesis B is a curved power law (3 free parameters): 

dN/dE = NoE-^" + '^^'''-^')\ (3) 

which corresponds to a parabolic law in a log(E'^ dN/dE) 
vs. logE representation. A parabolic shape is a nat- 
ural choice to describe the transition region of a break 
between two spectral components. By fitting the two 
functional forms (Eq. [5] and [3]) to the de-absorbed spec- 
trum one obtains values of the likelihood functions La 
and Lb- If hypothesis A is true the likelihood ratio 
R = — \n{LA/LB) is approximately distributed with 
one degree of freedom. One defines a probability 



P 



P (X^) dx^ 



(4) 



where p(x^) is the probability density function and 
Llmeas is the mcasurcd value of R. Hypothesis A will 
be rejected (and hypothesis B will be accepted) if P is 
greater than the confidence level, which is set to 95%. 

5. RESULTS 

We examined a wide range of redshifts values z be- 
tween 0.1 and 0.9 in steps of 0.02. Each time, the corre- 
sponding optical depth was calculated and the intrinsic 
spectrum of PG 1 553+113 was d etermined using the low 
limit model from iKneiske et al] ([2004). The probability 
of the likelihood ratio test and the intrinsic photon index 
Pint as a function of redshift z are shown in Fig. |4l 

The intrinsic photon index Fjnt as a function of the 
redshift is shown by the thick blue line in Fig. |4l A 2 cr 
confidence belt is drawn as blue shaded area. The re- 
sult without the normalization between the H.E.S.S. and 
MAGIC measured spectrum of PG 1553+113 is shown by 
the dashed blue line with a corresponding 2a confidence 
belt as grey shaded area. Assuming that the intrinsic 
photon index can not be harder than Fjnt = 1-5, we ob- 
tain a redshift limit of z < 0.69 with a confidence of 95%. 
Assuming the max i mally hard spectrum as proposed by 
iKatarzvhski et all ()2006l ) with Pint = 2/3, we obtain a 
redshift upper limit of z < 0.80. 

The intrinsic spectrum was considered to have a break 
if the likelihood ratio test gave more than 95% confi- 
dence. The probability of the likelihood ratio test as a 
function of the redshift is shown by the black points and 
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Fig. 4. — Constraints on the redshift of PG 1553+113 (see text 
for details). 

Black points and left Y-axis: probability of the likelihood ratio 
test (Eq.Hll. The black arrow indicates the minimum redshift {z = 
0.42), at which the break in the intrinsic spectrum of PG 1553+113 
is evident. 

Blue line and right Y-axis: intrinsic photon index Tint. Shaded 
areas correspond to 2 a confidence belt. The intrinsic spectrum is 
harder than this for z > 0.69, which is indicated by the blue arrow. 
The extreme case of P = 2/3 leads to an upper limit of 2 < 0.80 
as indicated by the dashed blue arrow. 

thick black line in Fig.[4l The smallest redshift, for which 
the test gave more than 95% confidence, is z = 0.42 (see 
Fig. [3]). The assumption that there is no break in the 
intrinsic spectrum of PG 1553+113 thus leads to an up- 
per limit on its redshift of 2; < 0.42. The result with- 
out the normalization between the H.E.S.S. and MAGIC 
measured spectrum is shown by the thin black line. The 
systematic uncertainty in the derived limits, arising from 
the applied normalization of 10% between the H.E.S.S. 
and MAGIC data, is below 10%. 

The test for the existence of a break cannot be ap- 
plied to the H.E.S.S. spectrum alone since the break oc- 
curs around the lower end of the energy region of the 
H.E.S.S. measurement. On the other hand, applying this 
test to the MAGIC data alone never reaches the required 
probability of 95% for the existence of the break due to 
insufficient statistics. 



6. CONCLUSION AND OUTLOOK 

In this study, we combined H.E.S.S. and MAGIC data 
using their good agreement and used a realistic minimum 
density of the EBL to reconstruct the intrinsic spectrum 
of PG 1553+113. We showed that the intrinsic photon 
index Pint becomes smaller than Pint = 1.5 at z = 0.69, 
which can be considered a robust upper limit on the red- 
shift of PG 1553+113. In case of the extreme emission 
scenario with Pint — 2/3, we obtain an upper limit of 
z < 0.80. Moreover, we showed that a break in the in- 
trinsic spectrum at about 200 GeV becomes evident at a 
redshift of z = 0.42. The break can either be interpreted 
as an upper limit on the redshift of PG 1553+113 or as 
evidence for a second emission component in the VHE 
spectrum of the object. The upper hmit of z < 0.42 
implies values of the intrinsic slope indicating that the 
peak of the high-energy component of the SED lies be- 
low few hundred GeV, as typically derived for the clos- 
est TeV blazars Mrk421 and Mrk501 (in low flux state). 
We note that increasing the statistics by combining the 
spectra of MAGIC and H.E.S.S. resulted in a moderate 
improvement of the redshift upper limit. On the other 
hand, the second method, which is based on the search 
for structures in the intrinsic spectrum and which re- 
sulted in much more stringent limits, became only feasi- 
ble using the combined data set. 

The knowledge of the distance of PG 1553+113 is cru- 
cial for the modeling of the emission processes in the ob- 
ject, especially for combining X-ray and VHE data. In 
case the redshift of PG 1553+113 is larger than z > 0.2 
it would not only be the farthest AGN detected in the 
VHE range so far but also by far the most energetic 
one. Further campaigns to determine the redshift of 
PG 1553+113 directly as well as extensive VHE obser- 
vations are encouraged. 
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